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A REPRESENTATION OF JACCHIA'S THERMOSPHERIC MODELS 
IN SPHERICAL HARMONICS 


P. Blum and I. Harris 


ABSTRACT 

The Jacchia models are represented in terms of spherical harmonic 
functions. This representation has the advantages of ease of comparison 
with theoretical and other observational models and data, mathematical 
analyticity and relative simplicity. The symmetry properties of the 
models are emphasized by this representation and some physical char- 
acteristics like the increase of the amplitude of the diurnal density 
variation with decreasing solar activity become more apparent. 
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A REPRESENTATION OF JACCHIA’S THERMOSPHERIC MODELS 


IN SPHERICAL HARMONICS 


The semi-empiric Jacchia thermospheric models (Jacchia 1965, 1970 and 1971) 
describe thermospheric densities as derived from satellite drag analysis. The 
densities are obtained in the models from a semi -empiric temperature profile 
and the global distribution of exospheric temperature. The exospheric tempera- 
tures are model parameters that are not necessarily identical with the true 
kinetic temperatures. 


In the Jacchia models the exospheric temperature depends upon the solar 
activity F, the average solar activity F during several solar rotation periods, 
the geomagnetic activity index K p , the day in the year, the local solar time 
(diurnal variation) and the latitude. The variation of exospheric temperature 
with local time and with latitude is given by a single mathematical expression; 
i.e., the distribution of exospheric temperature on the globe where the longitude 
is measured with respect to local noon. 


In Jacchia's model (1971) the exospheric temperature is given (in his notation) 
by: 


T c (1 + R sin m 6) 


( cos m 7? - sin m 6 

1 + R L 

\ 1 + R sin m 6 



(1) 


where T c depends only upon solar and geomagnetic activity. 

T c = 379° + 3° 24 F + l°3 (F - F) + 28°K p + 0°03 exp (K ) 

r = H + /3 + p sin (H + y) (- rr < r _< tt) 


1 



V 


\<t> - §1 

2 

\± + sj 

2 

where 0 is the latitude, S the declination of the Sun, H the LST counted from 
local noon, /3 = -37°, p - 6°, m = 2.2, R = 0.3, and y = 43°. 

This representation of the exospheric temperature includes absolute values; 
therefore it is not analytic. Especially its derivative does not exist at the poles 
(Blum and Harris, 1973). Jacchia's expression for the exospheric temperature 
distribution may be approximated by a sum of spherical harmonic functions. 

Such an approximation is useful because it allows a direct comparison with 
thermospheric models derived from other considerations like mass-spectrographic 
data (Hedin et al., 1973) or a theoretical treatment (Harris and Priester, 1962; 
Mayr, Harris and Volland, 1973). Furthermore, analytic expressions are handled 
with greater mathematical ease. The pressure gradients, which drive the thermo- 
spheric wind system, can only be calculated for the entire globe if the tempera- 
ture is an analytic function. 

In Jacchia's model the exospheric temperature is given as the product of the 
global minimum exospheric temperature, T c , and an expression for the tempera- 
ture dependence on latitude, local time and the declination of the sun; i.e. the 
day in the year. Thus 


T ~(F, F, K , 0, t A , tg) 


T c (F, F, K p ) • G (0, t A , t£ ) 


where tg is the local time and t A the day in the year. 


(2) 
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We have for the moment excluded the semi-annual variation as it will be treated 
separately. The above factorization of the exospheric temperature representation 
allows us to express the exospheric temperature in a simple, though very general, 
form as a sum of spherical harmonic functions P™ (x) where x is the cosine of 
the colatitude. 

E E “nm^A) P „"( X ) COS O t 0 + L(‘a) P n <» Sin ( m * E )) ( 3 ) 

For the numerical results of this paper we have used the usual normalization of 
the spherical harmonic functions (Jahnke-Emde, p. 111). is the mean global 
exospheric temperature that depends on F, F and K p . 

Expression (3) can be further simplified by making use of the following properties 
of Jacchia's model: 

(a) The phase of the diurnal variation is independent of latitude and season. 
Therefore, the cosine and sine terms of the various modes (i.e. the a 
and yS) can be combined so that all the P n m that have the same local time 
mode, m, have the same phase t . 

(b) The coefficients a mn and /3 mn depend only on the day in the year and can 
be expressed as a sum of Fourier modes with the argument the day in 
the year. These Fourier expressions become particularly simple if the 
days are counted from vernal equinox. Then only sine coefficients ap- 
pear in the symmetric terms (i.e. P™ with n - m = 2k) and cosine coef- 
ficients in the antisymmetric terms (n - m = 2k + 1). This property is 
due to the symmetry of the exospheric temperature distribution at equinox. 
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(c) The third harmonic of the diurnal variation has an amplitude that is 
less than 0.02 of the amplitude of the fundamental. This is less than 
the accuracy of drag data and has no physical significance in the Jacchia 
model. For this reason we have not included the third harmonic in the 
representation of the temperatures. Due to the non-linearity of the ex- 
ponential function the relative amplitude of the third harmonic of the 
density distribution is slightly enhanced. We have included it in the 
analysis of the density distribution although its physical significance 
for the densities is equally doubtful. 

4. A numerical analysis shows that no terms P™ with n > 4 occur with ampli- 
tudes larger than 5.10 -4 of the amplitude of the mean global term P Q °. Expres- 
sion (3) simplifies therefore to 

T ~ = E E ^(‘a) P n m ( x > cos m ( fc £ " O) ( 4 ) 

\ m-0 n =ra / 

with 

^(‘a) = E o (y n ( ^ s) cos k t A + y n ^ a > sink t A ^ (5) 

and y 0 ( 0 °' s > normalized to unity. 

The analysis also shows that no higher harmonics than the semi-annual appear 
in the Fourier analysis of the coefficients y nm (t A ). It should be noted that the 
semi-annual terms that appear in (5) arise from the variation of the declination 
of the sun and are quite small. These semi-annual terms should not be confused 
with the main semi-annual terms that are associated with the semi-annual effect. 
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Table 1 shows the numerical values of the various terms when the yearly average 
global term P Q ° is normalized to unity. Thus Table 1 is independent of solar 
activity. The relationship between the global minimum exospheric temperature 
T c in the Jacchia model and the global mean exospheric temperature used here 
is given by 

= 1- 1275 T c 

A synthesis of the exospheric temperatures based on Table 1 deviates less than 
1% from Jaachia's original representation. Considering the accuracy of drag 
data, no loss of physical significance results from the use of this representation 
as compared to the original model. 

THE SEMI-ANNUAL EFFECT AND THE GEOMAGNETIC ACTIVITY EFFECT 

In the Jacchia models the semi-annual and the geomagnetic activity correction to 
the exospheric temperature have no latitudinal dependence. They modify only 
the P 0 ° terms after the yearly average global term P 0 °, is normalized to unity. 

As the geomagnetic activity effect has no periodic time dependence, it is simply 
included in the global average T^ term and so causes no change in the numerical 
values of the coefficients given in Table 1. It should be noted that other empirical 
models of the thermosphere (Hedin et al., 1973) have deduced a latitudinal depen- 
dence of the geomagnetic activity effect. 

In the 1965 and 1970 Jacchia models the semi-annual variation was expressed as 
a correction to the exospheric temperature. We have Fourier analysed this cor- 
rection and obtained 
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( 6 ) 


T sa = F (0.2203 • cos (t A - 283.7 d) + 0.3309 • cos 2 (t A - 30 d) 

+ 0.0762 cos 3 (t A - 10.3 d)j 

The mean annual value of T sa is insignificant and has been neglected. The third 
annual harmonic has little physical significance, but we have included it in the 
semi-annual modification given in Table 1. The semi-annual temperature ampli- 
tude is directly proportional to the solar activity F, but the global mean exospheric 
temperature is only linearly dependent on F. To generalize the semi-annual cor- 
rection given in Table 1 for various values of solar activity we have to multiply 
the semi-annual amplitudes by a function of F. This function is approximated 
by F : 

J s a 


F 


s a 



- 0.124 



(?) 


In the Jacchia 1971 model the semi-annual variation is given in terms of densities 
and not temperatures; therefore it cannot immediately be included in a modifica- 
tion of Table 1. 

REPRESENTATION OF DENSITIES 

The thermospheric densities are obtained in the Jacchia models from the tem- 
peratures by adopting lower boundary densities and assuming diffusive equilibrium. 
Due to the non-linearity of the exponential function the relative amplitudes of the 
various spherical harmonic functions are altitude-dependent. There is also a 
dependence on solar activity. For this reason the direct representation of the 
densities as a sum of spherical harmonics is less useful then the representation 
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of the exospheric temperature because even the normalized coefficients depend 
on two parameters: altitude and solar activity. We have given the density repre- 
sentation at a solar activity of F = 150 for the altitudes 300, 400 and 500 km in 
Tables 2-4. These tables do not include the global semi-annual variation. 

Some spherical harmonic functions that did not appear in the analysis of the 
temperature, because their normalized amplitude was less than 5 • 10 ~ 4 (the 
threshhold we have chosen), will have larger amplitudes in the analysis of the 
densities due to the non-linearities and appear in the Tables 2-4. 

The decomposition of the global semi-annual density variation p gA (t A ) into 
Fourier modes according to the 1971 Jacchia model was given by Volland, Wulf- 
Mathies and Priester (1972). Our independent analysis gave almost identical 
results. Following Volland et al. the density correction at a height of 400 km is 
given by 

Ps\(*y) 

= 0.1034 cos (t y - 4d) + 0.1999 cos 2 (t y - 109 d) 

+ 0.0394 cos 3 (t y - 66 d) + 0.0188 cos 4(t y ~ 14 d) 

( 8 ) 

where t is the day in the year counted from January 1. 

Our results are 

PsA ( t s) 

— - 0.1048 cos (t A - 288 d) + 0.2025 cos 2 (t A - 28 d) 

+ 0.0399 cos 3 (t A - 107 d) + 0.0191 cos 4 (t A - 24.6 d) 

= 0.0255 cos t A - 0.1016 sin t A + 0.1155 cos 2 t A + 0.1664 sin 2 t A 

+ 0.0290 cos 3 t A - 0.0275 sin 3 t A - 0.0024 cos 4 t A + 0.0189 sin 4 t A 

(9) 
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Table 5 gives the density representation including the global semi-annual effect 
at a height of 400 km and a solar activity of F = 150. 

The representation of the drag data in terms of spherical harmonics given in 
this paper has the advantage of being analytic and relatively simple; i.e. using 
a small number of constants. It is transparent both from a mathematical and 
physical point of view and allows the evaluation of the physical significance of * 
some properties of drag data results. 
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Table 1 


Exospheric Temperature Distribution of the Jacchia Model 
(Independent of Solar Activity) 


(Excluding Global Semi-Annual Effect) 
Coefficients of Spherical Harmonics 


Spherical Harmonic 

Yearly Average 
cos (0) 

Annual 

cos t A sin t A 

Semi-Annual 
cos 2t A sin 2t A 

Third Annual 
cos 3t A sin 3t A 

P o° 

1.00 




P? 


.0528 



P 2 ° 

.0054 


-.0007 


P 3 


.0022 



P 4° 

.0041 


.0006 


p» 





Pj 1 COS (t £ - tj) 

.1300 


.0056 


P2 1 COS (t£ - t x ) 


.0022 



P 2 2 cos 2 (tg - t 2 ) 

.0077 





To include global semi-annual effect replace P 0 ° line by: 


p° 

0 

1.00 

.0056 F sa -.0326 F SA 

.0249 F sa .0427 F SA 

.0064 F sa -.0095 F sa 


t x = 14.4 
t 2 = 1.89 

t A is measured in days from the vernal equinox 
t £ is measured in hours from midnight 



Table 2 


Density Distribution of the Jacchia Model 
Height = 300 km F = 150 

(Excluding Global Semi-Annual Effect) 
Coefficients of Spherical Harmonics 



Spherical 

Harmonic 

Yearly Average 
cos (0) 

Annual 

cos t A sin t A 

Semi-Annual Third Annual 

cos 2t A sin 2t A cos 3t A sin 3t A 


Po° 

1.0 


-.0013 


Pi° 


.1416 

.0009 


P 2 ° 

.0126 


-.0025 


P 3 ° 


.0074 



P 4° 

.0112 


.0015 , 


po 




p/ 

cos (tg - tj) 

.3492 


.0150 

P 2 1 

cos (tg - tj) 


.0081 


P 3 1 

COS (tg - tj) 

.0013 


.0008 

p 2 

cos 2 (tg - t 2 ) 

.0217 


.0009 

P 2 

cos 2 (tg - 1 2 ) 




P 4 

cos 2 (tg - 1 2 ) 

.0012 



P 5 2 

cos 2 (tg - 1 2 ) 




P 3 3 

cos 3 (t - t 3 ) 

.0007 



*1 = 

14.40 

t 2 = 1.91 

t 3 = 3.5 



t A is measured in days from the vernal equinox tg is measured in hours from midnight 



Table 3 


Density Distribution of the Jacchia Model 
Height = 400 km F = 150 

(Excluding Global. Semi-Annual Effect) 
Coefficients of Spherical Harmonics 



Spherical 

Harmonic 

Yearly Average Annual 

cos (0) cos t A sin t A 

Semi-Annual Third Annual 

cos 2t A sin 2t A cos 3t A sin 3t A 


Po° 

1.000 


-.0018 


po 


.2091 

.0013 


p° 

2 

.0045 


-.0078 


P 3° 


.0104 

.0007 


P 4° 

.0163 


.0018 


po 




Pi 1 

COS (tg - tj) 

.5173 


.0219 

P 2 

COS (tg - tj) 


.0250 

.0007 

P 1 

COS (tg - t j) 

.0019 


.0007 

p2 

cos 2 (tg - t 2 ) 

.0396 


.0021 

P 3 

cos 2 (tg - t 2 ) 


.0010 


P 4 2 

cos 2 (tg - t 2 ) 

.0018 



P 5 2 

cos 2 (tg - t 2 ) 




P 3 3 

cos 3 (tg - t 3 ) 

.0009 



tl = 

14.39 t 2 

= 2.00 t 3 = 4.4 



t A is measured in days from the vernal equinox 


is measured in hours from midnight 



Table 4 


Density Distribution of the Jacchia Model 
Height = 500 km F = 150 

(Excluding Global Semi-Annual Effect) 
Coefficients of Spherical Harmonics 



Spherical 

Harmonic 

Yearly Average Annual 

cos (0) cos t A sin t A 

Semi-Annual Third Annual 

cos 2t A sin 2t A cos 3t A sin 3t A 


Po° 

1.000 

-.0022 


Pi 

.2726 

.0016 


P 2 ° 

-.0104 

-.0149 


P 3 ° 

.0121 

.0006 


P 4 ° 

.0211 

.0018 


P 5 ° 

.0017 


Pi 1 

COS (t g - t x ) 

.6759 

.0284 

P2 1 

COS (t g - t 3 ) 

.0476 

.0014 

p 1 
■*■3 

COS (tp - t x ) 

.0021 


p 2 

cos 2 (tp - t 2 ) 

.0604 

.0034 

P 2 

cos 2 (tg - t 2 ) 

.0025 


P4 

cos 2 (tp - t 2 ) 

.0024 


P 3 3 

cos 3 (tg - t 3 ) 

.0016 


tl = 

14.37 t 2 

= 2.05 t 3 = 5.1 


t A is 

measured in days from the vernal equinox tg 

is measured in hours from midnight 
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Table 5 


Density Distribution of the Jacchia Model 
Height = 400 km F = 150 

(Including Global Semi-Annual Effect) 
Coefficients of Spherical Harmonics 


Spherical 

Harmonic 

Yearly Average 
cos (0) 

Annual 

cos t A sin t A 

Semi- 
cos 2t A 

-Annual 
sin 2t A 

Third Annual 
cos 3t A sin 3t A 

P 0 ° 

1.000 

.0251 

-.1004 

.1122 

.1642 

.0286 

-.0271 

P° 

-.0106 

.0175 

.1973 

.0075 


-.0152 

.0135 

p0 

.0040 


-.0007 

-.0073 

.0007 



p° 


.0009 

.0099 



-.0008 

.0013 

P 4 ° 

.0164 

.0005 

-.0016 

.0036 

.0027 

.0005 

-.0005 

P 5 ° 








P* COS (tg - t x ) 

.5186 

.0136 

-.0511 

.0808 

.0851 

.0151 

-.0150 

P 2 COS (tg - t t ) 

-.0013 

.0021 

.0237 

.0009 


-.0018 

.0022 

P^ COS (tg - tj) 

.0019 



.0009 




P 2 2 cos 2 (tg - t 2 ) 

.0397 

.0010 

-.0039 

.0065 

.0065 

.0012 

-.0012 

P 3 2 cos 2 (tg - t 2 ) 



.0010 





P 4 2 cos 2 (tg - t 2 ) 

.0018 







P 3 3 cos 3 (tg - t 3 ) 

.0009 








t x = 14.40 t 2 - 2.00 t 3 = 4.4 

t A is measured in days from the vernal equinox tg is measured in hours from midnight 



